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ABSTRACT 

We examine electron-capture supernovae (ECSNe) as sources of elements heavier than iron in the 
solar system and in Galactic halo stars. Nucleosynthesis calculations are performed on the basis of 
thermodynamic histories of mass elements from a fully self-consistent, two-dimensional (2D) hydro- 
dynamic explosion model of an ECSN. We find that neutron-rich convective lumps with an electron 
fraction down to l^ im i n = 0.40, which are absent in the one-dimensional (ID) counterpart, allow for 
interesting production of elements between the iron group and N = 50 nuclei (from Zn to Zr, with 
little Ga) in nuclear (quasi-)equilibrium. Our models yield very good agreement with the Ge, Sr, Y, 
and Zr abundances of r-process deficient Galactic halo stars and constrain the occurrence of ECSNe 
to ~4% of all stellar core-collapse events. If tiny amounts of additional material with slightly lower 
y c ,min down to ^0. 30-0.35 were also ejected — which presently cannot be excluded because of the 
limitations of resolution and two-dimensionality of the model — , a weak r-process can yield elements 
beyond N — 50 up to Pd, Ag, and Cd as observed in the r-process deficient stars. 
Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abundances — stars: neu- 
tron — supernovae: general 



1. INTRODUCTION 

Electron-capture supernovae (ECSNe) have been pro- 
posed as possible origin of elements beyond iron, in 
particular of heavy r-process elements. This is due to 
the distinctive features of the collapsing O-Ne-Mg cores 
(|Nomotolll987h . The cores with their fairly small mass 
(~l.38.M0), steep surface density gradient, and dilute 
H/He envelope (in the super-asymptotic-giant-branch or 
SAGB stage) were expected to enable prompt explo- 
sions, where the ejection of neutron(n)-rich matter with 
an electron fraction (numb er of protons per ba ryon) of 
Y c « 0.14 could be possible (|Wanaio et al.ll2003| ). Chem- 
ical evolution studies also favor ECSNe at the lower end 
(~8-10M Q ) of the mass range of core-collapse super- 
novae (CCSNe) to acc ount for the abundance sig natures 
of Galactic halo stars (|Ishimaru k, Wanaiolll999D . 

Recent one-dimensional (ID) hydrodynamic simula- 
tions of O-Ne-Mg core collapse with elaborate neutrino 
transport , however, found delayed neutrino- driven ex- 
plosions ()Kitaura. Janka. fc Hillebrandtl 12 006). This re- 
sult was confirmed bv lDessart et al.l (120061) and ex cludes 
the prompt explosions of iHillebrandt eTaLl (fl98l . The 
shocked surface layers of the O-Ne-Mg core are also far 
from the r-process conditions envision ed by iNing et all 
(2007() as shown bv I Janka et"aT1 (|2008| ). 

Detai led nucleosynthesis studies bas e d on the ID re- 
sults of iKitaura. Janka. fc Hillebrandti ([2006), however, 
suggest that ECSNe can be the source of nuclei up 
to N — 50, in particular of Zn and light p-nuclei 
(up to 92 Mo), whose origins are not fully understood 
(Wanajo et a l. 200 9: see also Hoffman et al. 2008). 
iWanaio et alj (|2009D further speculated that the nucle- 
osynthetic outcome could be considerably altered in the 
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multi-dimensional case. 

Spectroscopic studies of Galactic halo stars indi- 
cate the presence of a significant fraction of r-process 
deficient Galactic halo stars that possess high Sr- 
Y-Zr (N = 50 ) abundances re l ative to heavy r- 
proce ss material (|Johnson fc Boltd 2002; H onda et al.l 
2006). This implies the contribution from an as- 
trophysical source making N = 50 species with- 
out heavy r-proces s elements dTravaglio et al.1 120041 : 
IWanaio & Ishimarul 120061: lOian fc WasserbursJ 120081). 
Recently I Arcones and Montesl ( 20101 ) . employing trajec- 
tories from ID hydrodynamic models and varying the 
corresponding Y" c , proposed that the baryonic wind blown 
off the surface of new-born ("proto-") neutron stars 
(PNSs) by neutrino heating could be production sites 
of such nuclei in the early Galaxy. 

New 2D explosion simulations of ECSNe exhibit n-rich 
lumps of matter being dredged up by convective over- 
turn from the outer layers of the PNS during the early 
stages of the explosion (§ 2; Muller, Janka, & Kitaura, in 
preparation), a feature that is absent in the ID situation. 
In this Letter we present the first nucleosynthesis study 
on the basis of such self-consistent 2D models. Without 
invoking any free parameters we show that ECSNe could 
be an important source of the elements between the iron 
peak and the N = 50 region in the Sun and the early 
Galaxy, perhaps even of light r-process nuclei up to Cd. 

2. TWO-DIMENSIONAL SUPERNOVA MODEL 

The nucleosynthesis analysis made use of about 2000 
representative tracer particles, by which the thermo- 
dynamic histories of ejecta chunks were followed in 
our 2D hydrodynamic calculation of an ECSN. The 
model was computed with a sophisticated (ray-by- 
ray-plus) treatment of the energy-dependent neutrino 
transport, using the Prometheus- Vertex code and 
the same microphysics (weak-interaction rates, nu- 
cle ar burning treatment, an d nuclear equation of state 
of lLattimer fc Swestvl Il991f) as in its ID counterpart 
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Fig. 1. — Snapshot of the convective region of the 2D simulation 
of an ECSN at 262 ms after core bounce with entropy per nucleon 
(s; left) and Y e (right). Mushroom-shaped lumps of low-Y e matter 
are ejected during the early phase of the explosion. 
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Fig. 3. — Eleme ntal mass fracti ons in the ECSN ejecta relative to 
their solar values (Loddcrs 2003), co mparing the 2D resu lts (red) 
with the ID counterpart (blue) from Wanaio ct al. (2009). Even- 
Z and odd-Z elements are denoted by circles and squares, respec- 
tively. The normalization band (see text) is marked in yellow. 
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Fig. 2. — Ejecta masses vs. Ye for the ID (blue) and 2D (red) 
explosion models. The width of a Y c -bin is chosen to be AY C = 
0.005. The minimum values of Y a are indicated for both cases. 



(jKitaura. Janka. "fc Hillcbrandt 2006). Some aspects of 
th e 2D model in comp arison to ID results were discussed 
bv lJanka et al.l (|2008f ). and a more complete description 
can be found in a separate publication (Miiller et al., in 
preparation) . 

The pre-collapse model of the O-Ne- Mg core emerge d 
from the evolution of an 8.8 Mq star (|Nomotol Il987[) . 
Because of the very steep density gradient near the core 
surface, the shock expands continuously, and a neutrino- 
powered explosion sets in at t ~ 100 ms p.b. in ID and 
2D essentially in t he same way and with a very similar 
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energy (~10 50 erg; Uanka et al 

In the multi-dimensional case, however, the negative 
entropy profile created by neutrino heating around the 
PNS leads to a short phase of convective overturn, in 
which accretion downflows deleptonize strongly, are neu- 
trino heated near the neutrinosphere, and rise again 
quickly, accelerated by buoyancy forces. Thus n-rich 
matter with modest entropies per nucleon (s ~13-15 fee; 
fee is Boltzmann's constant) gets ejected in mushroom- 
shaped structures typical of Rayleigh- Taylor instability. 
Figure Q] displays the situation 262 ms after bounce when 
the pattern is frozen in and self-similarly expanding. 

As a consequence, the mass distribution of the ejecta 
in the 2D model extends down to Y e m ; n as low as ~0.4, 



Fig. 4. — Mass fractions of the dominant isotopes of the even- 
Z elements vs. Y e for all investigated thermodynamic trajectories. 
Also shown are the mass fractions of 4 He (a-particles) (red) and 
the entropies of the trajectories (black dots, left vertical axis). 

which is significantly more n-rich than in the correspond- 
ing ID case (Y c ]^ in ~ 0.470. Figure [2] shows the Y c - 
histograms at the end of the simulations. The total 
ejecta masses are 1.39 x IO -2 Mq for the ID model and 
1.14 x 10~ 2 Mq in 2D, where the difference is partly 
due to the different simulation times, being ~800 ms and 
~400ms, respectively (core bounce occurs at ^50 ms). 
However, the ejecta after ^250 ms p.b. are only proton- 
rich, contributing merely to the Y e > 0.5 side in Fig. [2j 

3. NUCLEOSYNTHESIS FOR THE ECSN MODEL 

The nucleosynthetic yields are obtained with the reac- 
tion netw ork code (including n eutrino interactions) de- 
scribed in Wa naio et al.l (|2009f ) . Using thermodynamic 
trajectories directly from the 2D ECSN model, the cal- 
culations are started when the temperature decreases to 
9 x 10 9 K, assuming initially free protons and neutrons 
with mass fractions Y c and 1—1^, respectively. The 
final abundances for all isotopes are obtained by mass- 
integration over all 2000 marker particles. 

3 Note that the exact lower bound of the mass distribution vs. Y c 
in the ID case is highly sensitive to details of the neutrino trans- 
port, e.g. the number and interpolation of grid points in energy 
space. In a rec e nt sim ulation with improved spectral resolution, 
IHiidepohl eTaTI (5010) obtained Y ejmin = 0.487. 



Electron-capture supernovae as the origin of elements beyond iron 



3 



The res ulting elementa l mass fractions relative to so- 
lar values (jLod dcrs 2003]), or the production factors, are 
show n in Fig. [3] (red) com pared to the ID case (blue) 
from [Wanaio et al. (|2009l ). The "normalization band" 
between the maximum (367 for Sr) and a tenth of that is 
indicated in yellow with the medium marked by a dotted 
line. The total ejecta mass is taken to be the sum of the 
ejected mass from the core and the outer H/He-envelope 
(= 8.8M -1.38M +O.O114M = 7.43 M Q ). Note that 
the N = 50 species, 86 Kr, 87 Rb, 88 Sr, and 90 Zr, have the 
largest production factors for isotopes with values of 610, 
414, 442, and 564, respectively. 

As discussed bv IWanaio et al.l (j2009l ). in the ID case 
only Zn and Zr are on the normalization band, although 
some light p- nuclei (up to 92 Mo) can be sizably produced. 
In contrast, we find that all elements between Zn and 
Zr, except for Ga, fall into this band in the 2D case 
(Ge is marginal), although all others are almost equally 
produced in ID and 2D. This suggests ECSNe to be likely 
sources of Zn, Ge, As, Se, Br, Kr, Rb, Sr, Y, and Zr, 
in the Galaxy. Note that the origin of these elements 
is not fully understood, although Sr, Y, and Zr in the 
solar system are considered to be dominantly made by 
the s-process. The ejected masses of 56 Ni (— s> 56 Fe; 3.0 x 
10~ 3 M©) and all Fe (3.1 x 1 0~ 3 M ) are the sam e as in 
the ID case (2.5 x 10~ 3 M , IWanaTo et al.ll2009h . 

The fact that oxygen is absent in ECSN ejecta but 
a dominant product of more massive CC SNe, can pose 
a constraint on the frequency of ECSNe ()Wanaio et al.l 
2009). Considering the isotope 86 Kr with its largest pro- 
duction factor in our 2D model and assuming / to be the 
fraction of ECSNe relative to all CCSNe, one gets 



_J_ A ( 86 Kr)/Y ( 16 O) 
1-f Af( 8 6Kr)/M noEC (i6 ) 



0.050, 



(1) 



where A ( 86 Kr) = 2.4 x 10~ 8 and Y^( 16 Q) = 6.6 xlQ- 3 
are the mass fractions in the solar system ([Lodd crs 2003), 
M( 86 Kr) = 1.1 x lO~ 4 Af is our ejecta mass of 86 Kr, 
and M„oEc( 16 0) = 1.5 M the production of 16 by all 
other CCSNe, averaged over the stel lar initial mass func- 
tion between 13 M and 40 M (see IWanaio et alj|2009t 
INomoto et all 120061 ) . Equation (1) leads to / = 0.048. 
The frequency of ECSNe relative to all CCSNe is thus 
~4%, assuming that all 86 Kr in the solar system ex- 
cept for a possible co ntribution from the s-process (18%, 
lArlandini et~aI| [T999) , originates from ECSNe. This is in 
good agreement with the prediction from a recent syn- 
thetic model of SAGB stars (for solar metallicity models, 
IPoelarends et al.ll2008l) . 

The remarkable difference between the ID and 2D 
cases (Fig. [3]) can be understood by the element forma- 
tion in nuclear (quasi-)eqilibrium. Representative for all 
elements in the normalization band of Fig. [3l the final 
mass fractions of the isotopes 64 Zn, 74 Ge, 80 Se, 84 Kr, 
88 Sr, and 90 Zr (accounting for dominant fractions of 49%, 
36%, 50%, 57%, 82%, and 51%, respective ly, of their el- 
ements in the solar system, Loddcrs 2003) arc displayed 
for all tracer trajectories in Fig. |4j 

Nuclear quasi-equilibrium (QSE) makes nuclei heavie r 
than the Fe-group up to A ~ 90 (jMever et all 119981 ). 
64 Zn, 88 Sr, and 90 Zr are thus produced at Y e = 0.43- 
0.49, where the a-concentration (at the end of the calcu- 
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Fig. 5. — Elemental abundances for various Ye.min compared with 
the stellar abundance s of the r-process d eficient star HD 122563 
with [Fe/H] -2.7 JHonda et al.l |2006|) an d the r- process en- 
hanced star CS 22892-052 with [Fe/H] sa -3.1 HSneden et al.l2003T ). 
The abundances of HD 122563 and CS 22892-052 are vertically 



shifted to match Zr for Y r 



0.40 and Ba for Y. 



0.15, 



respecti vely. The Cd and Lu values in HD 122563 (open circles) 
arc from Rocdcrcr et al. (2010), and the Ge value (open square) is 
from lCowan et al.l 02005! ). 



lations) is Y( 4 He) = 0.001-0.1. The QSE with abundant 
a particles, however, is known to leave a deep trough in 
the abundance curve between A ~60 and 90 because of 
the strong binding at N = 28 and 50. This explains the 
substantial underproduction of elements around Z ^33- 
37 in the ID case (Fig. [3l blue line). 

As the a-concentration becomes small, QSE asymp- 
totes to nuclear statistical equilibrium (NSE) . Since NSE 
with neutron excess (Y e ~ 0.4) leads to nuclei heavie r 
than the Fe-group up to A « 84 ()Hartmann et al.|[l985|) . 
the trough can be filled by NSE-abundances assembled in 
the n-rich ejecta lumps. The small Y( 4 He) (< 10~ 5 ) at 
Y c ~ 0.40-0.42 (Fig. |4j is indicative of the consumption 
of almost all a-particles before freeze-out from NSE. Ac- 
cordingly, QSE asymptotes to NSE in this Y^,-range and 
yields substantial amounts of 74 Ge, 80 Se, and 84 Kr (made 
as 74 Zn, 80 Ge, and 84 Se), nuclei that cannot be created 
in q-ric h QSE. A similar re sult can be seen in the QSE 
study bv lMever et a l. ( 1998) with entropies and Y e -values 
close to ours here (see Fig. 14 in their paper). 

In the n-rich ejecta lumps NSE (or a-deficient QSE) 
conditions are established for several reasons. They have 
smaller entropies (s « 13-15 fee per baryon) than the 
other ejecta (where s ~ 15-20 fee per baryon; Figs. [TJ|4]). 
This favors a-particles to disappear when NSE ends 
as the temperature drops. In addition, the a's be- 
come easily locked up and tightly bound in nuclei, i.e., 
their separation energies are large (cf., e.g., Fig. lb in 
iWoosley Sz Hoffmanlll992l ). because nuclei with n-excess 
do not readily release a's to move farther away from j3- 
stability. 

We find no sign of r-processing in the n-rich lumps. 
Our present calculations are limited to the first <400 ms 
after bounce and do not include the neutrino-driven PNS 
wind. The latter, however, turned out to have proton ex- 
cess in ID models of t he long-term evolution of ECSNe 
(|Hudepohl et al.lfeOlOl) . It thus makes only p-rich iso- 
topes with small production factors (below unity, Wanajo 
et al., in prep.) and has no effect on the discussed results. 

4. ECSN YIELDS AND GALACTIC HALO STARS 
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Our results in § 3 also imply that ECSNe can be the 
source of Sr, Y, and Zr as observed in r-process deficient 
Galactic halo stars (Fig. [5]). A number of such stars with 
detailed abundance determinations indicate, however, a 
possible link with the elem ents beyond N = 50, e.g. Pd 
and Ag (|Honda et al.ll2006t l. Our ECSN models cannot 
account for the production of such elements, but in their 
ejecta a small c hange of Y e can dr astically change the 
nucleosynthesis (|Wanaio et a l. 2009). Due to limitations 
of the numerical resolution and the lack of the third di- 
mension, or some sensitivity to the nuclear equation of 
state, it cannot be excluded that ECSNe also eject tiny 
amounts of matter with Y" c ,mm slightly lower than pre- 
dicted by the 2D simulation. 

We therefore compare the nucleosynthesis for Y Ciin i n = 
0.40 of our ECSN model and for artificially reduced 
values of F e>mill = 0.35,0.30,0.25,0.20, and 0.15 with 
the abundance patte rns of representati ve r-process de- 
ficient (HD 12 2563. iHonda et all [2006ft and enhanced 
fCS 22892-052. ISneden et all 120031) stars (Fig. 1ST) . For 
that we use the thermodynamic trajectory of the lowest 
Y e (= 0.404) of the original model but apply Y c down 
to 0.15 in steps of AY C = 0.005. The ejecta masses in 
these additional Y^-bins are chosen to be constant with 
AM = 2 x 10" 5 M Q in the cases Y" 0imin = 0.35 and 0.30, 
and AM = 10" 5 M Q for the other Y 0!min . 

Figure [5] shows that Y e < 0.35 is needed to obtain el- 
ements beyond N — 50. A remarkable agreement with 
the abundance pattern in HD 122563 up to Cd (Z = 48) 
can be seen for Y e ,min = 0.30. Such a mild reduction of 
l^min in the ECSN ejecta is well possible for the reasons 
mentioned above. A reasonable match of the heavier part 
beyond Z = 48 requires i^, m in ~ 0.20. This, however, 
leads to a poor agreement for Ag and Cd. We there- 
fore speculate that ECSNe could be the sources of the 
elements up to Cd in r-process deficient stars, and the 
heavier elements are from a different origin. Moreover, 
^e,min = 0.15 is necessary to reproduce the abundance 
pattern of r-process enhanced stars like CS 22892-052. 
Such a low Y seems out of reach and disfavors ECSNe 
as production sites of heavy r-process nuclei. 

We note that the neutron-capture reactions start 
from seeds with A ~ 80 formed in nuclear equilib- 
rium (§ 3), not by an a-capture process. We therefore 



prefer to call the described process producing the ele- 
ments be yond N = 50, presumably up to Cd, "weak r- 
process (|Wanaio fc Ishimarul 12009: IHonda et al.ll200l 
rather than q-process or charged-particle proc ess 
(IWooslev fc Hoffman! Il992 IQian fc Wasserburdl2008h . 

5. SUMMARY AND OUTLOOK 

Using ejecta-mass tracers from a self-consistent 2D ex- 
plosion model we computed the nucleosynthesis of EC- 
SNe. Because of a brief phase of convective overturn and 
the very fast explosion that is characteristic of O-Ne- 
Mg cores, n-rich lumps with Y e down to 0.4 are ejected. 
These allow for a sizable production of the elements from 
Zn to Zr in nuclear (quasi-)equilibrium. The model yields 
Ge, Sr, Y, and Zr in very good agreement with abun- 
dances of r-process deficient Galactic halo stars. A mild 
reduction of the minimum Y c to ^0.30-0.35, which can- 
not be excluded due to limited numerical resolution and 
the lack of the third dimension, leads to a weak r-process 
up to the silver region (Pd, Ag, and Cd), again well 
matching these elements in r-process deficient stars. The 
formation of heavy r-process nuclei requires Y c to be as 
low as ~0. 15-0.20 and seems out of reach for our models. 

We therefore determine ECSNe as an important source 
of Zn, Ge, As, Se, Br, Kr, Rb, Sr, Y, and Zr in the 
solar system and the early Galaxy. Our models, however, 
underproduce Ga, which suggests the s-process as the 
origin of this element. The frequency of ECSNe is thus 
constrained to ~4% of all CCSN events on average over 
the Galactic history, but could have been higher at early 
Galactic epochs, compatible with the commonality of r- 
process deficient halo stars. 

Future better resolved and in particular 3D models will 
have to elucidate the role of ECSNe as site of the weak 
r-process. Also important are new abundance studies of 
r-process deficient stars for more complete information 
on the elements from Zn to Zr and on weak r-process 
products between Nb and Cd. 
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